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We use an almost model-independent analytical parameterization for pp and pp elastic scattering 
data to analyze the eikonal, profile, and inelastic overlap functions in the impact parameter space. 
Error propagation in the fit parameters allows estimations of uncertainty regions, improving the 
geometrical description of the hadron-hadron interaction. Several predictions are shown and, in 
particular, the prediction for pp inelastic overlap function at ^/s — 14 TeV shows the saturation of 
the Froissart-Martin bound at LHC energies. 



I. INTRODUCTION 

o 

T-H The general belief that Quantum Chromodynamics, QCD, is the correct theory of strong interactions being the 

^^ only one able to describe this physical regime. However, reliable predictions can actually be extracted for a very 
^^ restricted class of physical quantities, i.e., in general QCD can only handle hard processes (perturbative processes). 
>^ Soft processes (high-energy elastic hadron scattering, for example) are usually described by models and some of them 
^ are "QCD inspired" |il|. However, there are models reaching remarkable results as, for example, Cheng and Wu model 
^> [2] (pre-QCD) that asymptotically predicts the total cross section, atot{s), will increase as R{s) where R{s) is a typical 

interaction radius which grows logarithmically with energy, atotis) ^ log s. 

^"^ Actually, explain high-energies pp and pp elastic scattering is a very difficult task since it is necessary to understand 

the formally theoretical results obtained within Axiomatic Quantum Field Theory, as well as the experimental available 

I— I data from hadron-hadron collision experiments realized over decades. In this way, models of elastic scattering must 

i-G obey theorems and bounds as the Froissart-Martin bound, Pumplin bound, Pomeranchuk theorem, and others since 

»-^ these theoretical results are model- independent. On the other hand, these formal results do not allow us to write a 

O functional form for pp and pp elastic scattering amplitude. Consequently, the experimental data play a central role in 

QJ this kind of analysis. Nowadays, experimental data about pp and pp elastic scattering allows a novel understanding 

i-G of scattering reactions because its data set are the most complete and that reaches the highest energies at colliders. 

The total cross section is an example of the importance of experimental data analysis since before the advent of ISR 

^— I collider, it was assumed that the pp and pp cross sections were approaching a common constant value of about 40 mb, 

J> at y/s — >■ oo. When the ISR was turned on, some results showed that pp cross section was rising with energy, more 

r ■ precisely, had a rise that could be fitted with a log s leading term. 

t ■ In [3l a phenomenological approach to explain high energies pp and pp total cross sections, differential cross sections 

and p parameter is presented. This parameterization is based on bounds and theorems of the Axiomatic Quantum Field 
Theory and on pp and pp experimental available data above -y/s — 20 GeV for the total cross sections, differential cross 
^O sections, and p parameter. This almost model-independent approach may be improved with additional theoretical 
information as well as using future experimental measurements at RHIC, LHC, and TEVATRON colliders allowing 
more precise descriptions. Furthermore, predictions in the momentum transfer space were performed and compared 
with those in the literature. In [S] a more detailed discussion is presented. 

On the other hand, there exist a classical description of the scattering amplitude based on the impact parameter 
representation, which the corresponding amplitude is perpendicular to the beam direction and thus is the same in the 
5_^ laboratory and center-of-mass system. The elastic scattering amplitude in the momentum transfer space allows, in 

d the impact parameter space, an estimate of the interacting particle sizes and an analysis of their internal structures. 
There are several papers about this theoretical approach in the literature with pros and cons [4j. In particular, Amaldi 
and Schubert had performed a comprehensive and coherent analysis oi pp data at ISR energies, which had shown that 
the rise of the total cross section is a peripheral phenomenon |S] . 

Using the well established parameterization taken from [J we investigate here the eikonal, profile, and inelastic 
overlap functions using the impact parameter representation via the Fourier-Bessel transform of the elastic scattering 



in 



o 
o 



X 



Permanent Address. 



^Electronic address: sergiodc@ufscar.br.br 



Phenomenological Analysis of pp and pp Elastic Scattering 2 

amplitude. Numerical values of the fit parameters used here were obtained in the momentum transfer space from 
simultaneous fits for the total, differential cross section and p parameter and then used to describe eikonal, profile and 
overlap functions in the impact parameter space. Standard error propagation in the fit parameters is also performed 
allowing more realistic descriptions of these functions. 

This paper is organized as follows. In section [TT} the basic formalism of the scattering amplitude in the impact 
parameter representation is presented. In section |lll| we show the numerical values of the parameters used here and 
present the results for the eikonal, profile and inelastic overlap functions. In Section [TV] predictions to other energies 
not used in the fitting procedures are shown, and in section IV] we present our final remarks. 

II. BASIC FORMALISM 

We consider elastic scattering amplitude in a simple geometrical picture using an impact parameter b and the 
Fourier-Bessel transform of the transferred momentum ~t — q^. In this picture, the profile function is introduced in 
high energy diffractive processes to describe the shape of the collision partners in the plane transverse to the beam 
axis. Assuming that longitudinal momentum transfer is negligible the profile function can be written in the impact 
parameter representation as [3 [TT] 



F{s, q')^i I bdbJo{qb)[l - e'^^''''^] = i I bdbjQ{qb)T{s, b) 



(1) 



where \/s is the energy in the center-of-mass system, Jq is a zeroth order Bessel function, and x(s, b) is the eikonal 
function written as 



x{s, b) — Rex(s, b) + ilmx(s, b). 



(2) 



The eikonal correspond to the continuum (complex) phase shift and in both the high energy limit and the semi- 
classical approximation, x(s, b) — 25{s, b). Moreover, Imx(s, b) is the so-called opacity function related to the matter 
distribution inside the incident particles. Explicitly, the real and imaginary parts of the eikonal function can be 
expressed in terms of the real and imaginary parts of the profile function r(s, b) as 



Rex(s, b) — tan 



Imr(s, b) 
_Rer(s,6)-l 



(3) 



Imx(s, b) = In 



Wr(s,6) + [1-Rcr(s,6)]2 



(4) 



The unitarity condition in impact parameter representation sets a link among total {(Jtot), elastic (cei) and inelastic 
(o'inei) differential cross sections. Consequently, it gives a remarkable opportunity to study at the same time both 
elastic and inelastic properties of the interaction. These physical quantities can be written in terms of r(s, b) as |11) 



ae/(s)= dH\T{s,b)\' 



(5) 



a.neiis) = f d^b[2Rer{s,b) - \r{s,b)\% 



(6) 



cTtotis) = 2 dHReT{s,b). 
In the impact parameter representation unitarity condition can be written as 



(7) 
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2Rcr(s, b) = \T{s, 6)p + G,„e/(s, b), (8) 

where Ginei{s,h) is the (real) inelastic overlap function. We can relate both overlap and inelastic cross section by 

/■OO 

0-me;(s) = / (fbGinel{s,b). (9) 

Jo 

In (Is]), |r(s, 6)p represents the shadow contribution of elastic channel, and it connects both elastic and inelastic 
channels. Using the constraint linx{s, b) > 0, it is easy to see that Ginei{s, b) represents the probability of an elastic 
scattering in (s, b) given by 



G,nelis,b) = l-e-^''-^^^''>^<l. 



(10) 



In general, derivation of the above formulas follows the usual argument which tends to imply that the impact 
parameter is a high energy, small angle approximation. However, using the Watson-Sommerfeld transform, Islam 
proves that the impact parameter representation is an exact representation of the scattering amplitude, valid for all 
energies and angles [12] • 

The basic assumptions about the almost model-independent analytical parameterization for pp and pp elastic 
scattering amplitude was well established in [3 . Besides those assumptions, we assume that F{s, q^) is also valid to 
g^ — )■ oo obtaining a unique r(s,b) (Appendix A). 

In the momentum transfer space, the elastic scattering amplitude for pp and pp can be written as 



s 






^21 ^-0.(sW 



ai{s)e 



ai{s)e 



^ + i'^a{s)t 



^/3(s)«' 



(11) 



i=l 



s 



1=1 



ai(s)e 



aj(s)e 



I +i^a{s)i 



-Ks)q^ 



(12) 



i=l 



where K is the subtraction constant, a^ (s) = A^ + Bi In s + Q In^ s, a^ (s) = Ai + i?^ In s + Ci In^ s, Pi (s) = D-i + S^ In s, 
Pi{s) = Di + Ei Ins, and Ai, Ai, Bi, Bi, C'i, Gt, Di, Di, Ei, and Ei are fit parameters. 

Using the Fourier-Bessel transform, the above elastic scattering amplitude provides the real and imaginary parts of 
the pp and pp profile function 



Rer(s, b)pj, = j-Y. 



ai{s) 



(13) 



Imr(s, b)^ 



a^_a^(sM^ 
8^L2/3.(s) 8/?3(s) 



E 



(4A(s) - b^ 



e43.(a 



32 



E 



ai{s) 
2A(s)' 



,43711) 



ai{s) 
2ft (s)' 



,43i(s) 



(14) 



Phenomenological Analysis of pp and pp Elastic Scattering 



Rer(.,.),,= -i-g|W,Sfe, ,15) 



1 ^ 

ImT{s,b)pp ^--^Yl 



«K^) ^^(^)^'^(^\^barP^^s)-b') 



3 r 



eVi(s) + 



--E 

2—1 



'il—Lpifiiis) Il—Lgnii(s) 

2ft (s) 2ft (s) 



(16) 



Notice that the real (imaginary) part of the profile function is connected with the imaginary (real) part of the elastic 
scattering amplitude. Therefore, if at high-energies, s — >■ oo, the elastic scattering amplitude (momentum transfer 
space) is imaginary dominant, then the profile function (impact parameter space) will be real dominant. 

In the impact parameter representation, the blackness growth speed is about zero at 6 = 0, leading saturation of 
the unitarity limit. Therefore, if imaginary part of the profile function can be neglected, then imaginary part of the 
scattering amplitude tends asymptotically to a limit, i.e., it cannot grow arbitrarily fast with the energy increase. 
As well established, the optical theorem in the momentum transfer space connects the imaginary part of the forward 
elastic scattering and total cross section, hence it is possible at future LHC energies obtain clear indications of 
saturation of the Froissart-Martin bound. 

III. FIT RESULTS 

The efficient description of the global clastic scattering data with an economical number of free parameters is a 
very difficult task since almost relevant physical quantities depend on both real and imaginary parts of F{s,q'^). 
Furthermore, model free parameters are, in general, so intricate that practically avoid identifying the explicit role 
of each parameter in the description of the experimental data and, in most of the cases, numerical calculations may 
introduce bias or loss of information, and does not allow standard error propagation. Therefore, the use of empirical 
or model-independent information about the physical quantities that characterize the elastic scattering amplitude is 
crucial to minimize these problems. 

Keeping in mind the above explanation, the almost-model independent parameterization proposed in [3] reaches 
global and efficient description of the elastic scattering data. This parameterization, although entail many parameters 
(30 free parameters to be exact), allows a global analysis of all elastic experimental data {atot, P, and differential cross 
section), allowing us to determine the uncertainties of the fit parameters and take into account the error propagation 
with standard methods. Furthermore, simultaneous fits for pp and pp experimental data in the forward direction 
as well as including data beyond the forward direction may result in a predictive formalism in both energy and 
momentum transfer f3]. 

The numerical values of the fit parameters obtained in the momentum transfer space are used to describe the profile, 
eikonal, and inelastic overlap functions in the impact parameter space. It is expected that the predictive character of 
these parameters remain in the impact parameter representation. In Table [l] the fit parameters obtained from data up 
to qfnax = 2.0 GeV^ and up to q'^^^x — 14.0 GeV^ are displayed. For the sake of the simplicity, only profile, eikonal, 
and inelastic overlap functions achieved from data up to q^ax = 14.0 GeV^ are shown. 

If we use the Fourier-Bessel transform, then the elastic scattering amplitude in the momentum transfer space implies 
in analytic functions in the impact parameter space. Therefore, error propagation in the fit parameters allows the 
estimation of the uncertainty regions [14], i.e., the uncertainties in the free parameters have been propagated to the 
profile, eikonal and inelastic overlap functions in the impact parameter space. We may estimate the confidence region 
associated with all these functions by adding and subtracting the corresponding uncertainties. We choose not to 
show the error bands as thatched areas since there are sometimes very narrow bands indeed turning this kind of 
representation useless here. 

A. Profile function 

The real and imaginary parts of the profile function for pp and pp are shown in Figs. [T] [2] and [3] The real part 
of r(s, b) is dominant, thus elastic scattering amplitude is almost imaginary at large momentum transfer. Therefore, 
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TABLE I: Results to pp and pp simultaneous fits to da/dq , atot, p at q^ax = 2.0 GeV 
K = -0.1053 ± 0.0048 and K = 49.7 ± 1.7, respectively. All parameters in GeV~^ K 
Ci + C2 + Ca — C2 — Ca • 



and qmax ~ 14.0 GeV , where 
= -0.1053 ± 0.0048 and Ci = 





pp 






pp 




ql,, = 2.0 GeV2 


ql,, = 14.0 GeV2 




ql,, = 2.0 GeV2 


ql,, = 14.0 GeV2 


Ai 


91.13±0.28 


109.70±0.28 


Ai 


119.61±0.44 


112.28 ±0.44 


Bi 


-12.939±0.039 


-16.529±0.039 


Bi 


-2.486±0.073 


-0.468 ±0.074 


Ci 






Ci 


-0.0174±0.0038 


-0.1673 ±0.0039 


Di 


-7.79±0.33 


-8.91±0.32 


Di 


3.134±0.067 


3.170 ±0.069 


El 


2.908±0.051 


3.045±0.050 


El 


0.4884±0.0078 


0.4860 ±0.0082 


A2 


16.82±0.22 


-4.06±0.23 


A2 


14.51±0.15 


10.23 ±0.15 


B2 


7.071±0.030 


11.387±0.030 


B2 


-6.730±0.024 


-6.756 ±0.027 


C2 


0.3027±0.0047 


0.0952±0.0047 


C2 


0.9035±0.0027 


0.9613 ±0.0029 


D2 


1.647±0.014 


1.290±0.014 


D2 


-1.549±0.011 


-1.476 ±0.013 


E2 


0.4646±0.0022 


0.5097±0.0023 


E2 


0.5521±0.0011 


0.5645 ±0.0012 


^3 


0.2582±0.0049 


1.0554±0.0077 


A3 


-17.160±0.055 


-8.148 ±0.031 


i33 


-0.09894±0.00074 


-0.3607±0.0013 


B3 


2.6083±0.0060 


1.2313 ±0.0040 


C3 


0.5921E-02±0.0070E-02 


0.02372±0.00012 


C3 


-0.11298±0.00038 


-0.05572 ±0.00025 


^3 


0.303±0.0019 


0.6454±0.0081 


D3 


1.5672 ±0.0084 


0.9272 ±0.0072 


i?3 





0.0176±0.0011 


E3 





0.03859 ±0.0011 



the elastic scattering is mainly the shadow of the inelastic channels and the momentum transfer dependence reflects 
the shape of the interacting hadrons [13] . 



Eikonal function 



The imaginary part of the eikonal, the so-called opacity, can be viewed in Fig. |4] The real part of the eikonal is 
always positive, which means attenuation and can only come from particle production diagrams llj not considered 
in our approach. 

It is interesting to note while not the purpose of this work that the eikonal in the impact parameter space, x(s, b) 
can be connected to eikonal in the momentum transfer space, x{s,q'^) by 



X{s,q )= / bdbJo{qb)x{s,b), 
Jo 

then x(s,5^) allows possible connections with quantum field theory since elementary cross sections are expressed in 
the momentum transfer space as well as form factors of the nucleons [S]. In general, the inputs for x(s, q^) are based 
on analogies with geometry, optics, or microscopic concepts in well established QCD bases. Of course, the crucial test 
for any model comes from the experimental data on the physical quantities that characterize the elastic scattering. 



C. Inelastic overlap function 



Amaldi and Schubert pointed out [5] that a precise knowledge of the real part of the scattering amplitude is not 
required for the impact parameter analysis since the effect of the real part on Ginei{s^b) is negligible. The inelastic 
overlap function can be viewed in Figs. [5] |6] and [7J In Fig. |8] are displayed the inelastic overlap functions to pp at 
y/s = 52.8 GeV, ^s = 7.0 TeV, y/s = 14.0 TeV, and to pp at ^ = 546.0 GeV, and at ^i = 1.8 TeV. Notice that 
the overlap function grows with the energy increase, but it is a non-increasing function of b for all smaller values of 
b. The Dubna Dynamical Model [S] of hadron-hadron scattering at high energy is based on the structure of a hadron 
as a compound system with a central region which the valence quarks are concentrated and log-distance region filled 
with a color-singlet quark-gluon field. As a result, the hadron amplitude can be represented as a sum of a central and 
a peripheral part and predicts the saturation of the overlapping function starting at ^/s ~ 1.8 TeV, as obtained here. 
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Contrary to Kundrat and Lokajicek work [51 [TU], the inelastic overlap showed here will be non- negative for any 



value of b. That occurs because the scattering amplitudes (11) and (12 1 are non-oscillating functions around zero at 
high values of b. Furthermore, Kundrat and Lokajicek define their inelastic overlap function, GfJ^i{s, b), using unitary 
condition in terms of the impact parameter amplitude 

, / M V{s, b) exp(iRex(s, b)) - 1 
hei(s,b) = — , 

where hei{s,b) ~ exp(— Imx(s, 6)). On the other hand, here Gineiis,b) is defined using unitary condition in terms of 
the profile function 

r(s, 6) = 1 - ?/(s, b) exp(zRex(s, b)) 

resulting Ginei{s,b)^^ — Ginei{s,b)/4:. The saturation of the Froissart-Martin bound may not be reached at LHC 
energies using this approach. 

IV. PREDICTIONS 

Figs, ^and 10 show predictions for |r(s, 6)p, 2Rer(s, b), and Ginei{s, b) at energies not available in the fit procedure 
described in 3 . 

At ^/s = 14 TcV, the inelastic overlap function almost reaches the saturation, i.e., Ginei{s,b) « 1. From (10 1, 
we may conclude that Imx(s, b) ^ oo and so the total cross section indicates the saturation of the Froissart-Martin 
bound. This saturation at y's — 14.0 TeV may imply corrections in the formalism used to describe this kind of 
scattering process. Therefore, some novel physical mechanism must be used to establish the correct growth of the 
elastic total cross section, if LHC run confirms this prediction 



On the other hand, |r(s, 5)|2 = [Rer(s, b)]'^ + [lmT{s, 6)]^ and from Fig. [lO]we see that Rer(s, 6) w 1 and Imr(s, b) 
is very small at 6 = 0. Therefore, the profile function is real dominant, i.e., the elastic scattering amplitude in the 
momentum transfer space will be imaginary dominant. Consequently, we expect a negligible real part of the scattering 
amplitude at ^/s = 14 TeV for q^ small values. 

V. FINAL REMARKS 

Empirical parameterization of the elastic scattering amplitude and fits of the available data have widely been used as 
a source of model-independent analysis and constitute important tactics that can contribute with the establishment of 
novel theoretical calculational schemes. The approach here applied is used to determine several quantities of interest, 
such as the profile, the eikonal, the inelastic overlap functions. With some additional hypotheses, even information 
about the form factors (momentum transfer space) can be extracted. 

Results in the impact parameter representation were obtained using the Fourier-Bessel transform of the elastic 
scattering amplitude in the momentum transfer space. In the impact parameter representation, we obtain the profile, 
eikonal, and the inelastic overlap functions in a confidence interval due to standard error propagation in the fit 
parameters. All results in the parameter space representation show that the maximum of their distributions lies at 
positive values of the impact parameter b. In the parameterization proposed in [3] an important assumption is that 
real part of the scattering amplitude may be neglected in comparison to the imaginary part at sufficient high values 
of transferred momentum. This assumption implies in a peripheral character to the pp and pp elastic processes, i.e., 
it occurs at small b (central character). 

In the almost model-independent analysis performed here, Ginei{s,b) is practically saturated at y/s — 14 TeV. 
Therefore, future experimental data obtained at LHC could be compared with this result allowing a confidence test 
of the parameterization. Furthermore, this premature saturation may show some novel physical mechanism in this 
energy region. 

As a remark, we note that the impact parameter description does not claim to be able to fit the data quantita- 
tively. However, it can provide evidence of fundamental features of high-energy scattering contributing for further 
developments of more elaborate diffractive models. 
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Appendix A: Uniqueness of the Profile Function 

Consider the known elastic scattering amplitude F{s, q^) written as a Fourier-Bessel transform in the range (0 < 

/"OO 

2 



F{s,q^)^i bdbjQ{qb)T{s,b). 
Jo 

The profile function is now known. On the other hand, consider the unknown amplitude F(s, q^) with respective 
unknown profile function 

/>oo 
r{s,b)=t q'dq'Uhq')F{s,q'^). 



The combined profile function is r(s, h) -\- r(s, h) and for (0 < g^ < q^ax) ^^ obtain the known amphtude 

POO 

F{s,q^) = i bdbJoiqb) [r{s, b) + f(s, b)] . 
Jo 

On the other hand, taking {q^ax < q^ < oo) we obtain 

/■OO 

F(s, q^) = l bdbJo{qb)[r{s, b) + f(s, b)l 



the unknown amplitude. Of course, we may construct an arbitrary number of scattering amplitudes (or profile 
functions) each one well established in its range. Therefore, for the sake of simplicity we assume F(s, q^) = F{s, q^) 
and hence, the profile function obtained in this work is unique. 
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1: Real e imaginary parts of the profile function to pp at ^/s 
J. — 14.0 GeV^. Dashed lines represent the error propagation. 
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FIG. 3: Real e imaginary parts of the profile function to pp at ^/s = 1800.0 GeV. Results obtained from the fit parameters to 
14.0 GeV^. Dashed lines represent the error propagation. 
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FIG. 4: Imaginary part of the eikonal function to pp at ^/s — 52.8 GeV (left) and to pp at ^/s = 546.0 GeV (riglit). Results 



obtained from the fit parameters to q 



14.0 GeV . Dashed lines represent the error propagation. 
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FIG. 5: Inelastic overlap function to pp at ^/s — 52.8 GeV (left) and at ^/s — 62.5 (right). Results obtained from the fit 



parameters to q 



14.0 GeV . Dashed lines represent the error propagation. 



Phenomenological Analysis of pp and pp Elastic Scatterin, 

1 



13 



^inel 



Gi 



nel 




FIG. 6: Inelastic overlap function to pp at ^ = 53.0 GeV at q^^^^ = 14.0 GeV^ (left) and at ^ = 546.0 GeV (right). Dashed 
lines represent the error propagation. 
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FIG. 7: Inelastic overlap function to pp at ^/s — 1.8 TeV to q'^rnax ~ 14.0 GeV^. Dashed lines represent the error propagation. 



0.4 



--_ .-^ 




- 








\ ■■''V 






\ ■ ■■■■■>^ 




■ 




^ 


■ 


\ '■ '^ 


■^c. 


- 


; \ 




\ 


■ 






■ 




\---\/\ 


- 





I i r^^-7~r ~r 1 


"^-^-^ 



1.5 
b(fm) 



2.5 



FIG. 8: Inelastic overlap function to: pp at ^ = 52.8 GeV (solid line) (ISR), pp at ^s = 546.0 GeV (dotted line) (SPS 
Collider), pp at ^ = 1.8 TeV (dashed line) (Tevatron), pp at ^i = 7.0 TeV (dash-dotted line) (LHC), and pp at ^ = 14.0 
TeV (long dashed line) (LHC). The inelastic overlap function grows with increasing energy at each b. 
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FIG. 9: Comparison between Ginei{s,b) (dotted line), 2Rer(s, 6) (dashed line) and |r(s,6)|^ (solid line) to pp at ^/s =200.0 
GeV. Results obtained from the fit parameters to q^^ax = 14.0 GeV^. 




FIG. 10: Comparison between Gi„ei{s,b) (dotted line), 2Rer(s,&) (dashed line) and |r(s,6)|^ (solid line) to pp at ^/s =14.0 
TeV. Results obtained from the fit parameters to q^ ^nax ~ 14.0 GeV^. 



